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Multi-constraints 


When  there  are  time-varying  parameters  in  motion  equations,  the  optimal  guidance  law  with  multi¬ 
constraints  generally  cannot  be  solved  analytically.  Based  on  block  pulse  functions,  a  design  method  of 
optimal  guidance  law  is  presented  combining  optimal  control  theory  and  numerical  value  computation 
for  time-varying  systems.  The  presented  guidance  law  can  optimize  the  combination  of  landing  angle, 
miss  distance,  and  control  energy  consumption.  Using  both  the  proposed  guidance  law  and  the 
proportional  navigation  law,  ballistic  simulations  are  made.  Compared  to  the  proportional  navigation  law, 
the  optimal  guidance  law  is  able  to  more  than  double  the  landing  angle.  Because  of  the  steep  terminal 
trajectory,  the  strike  accuracy  and  damage  effects  are  increased.  Averaging  the  time-varying  coefficients 
of  the  optimal  guidance  law,  a  suboptimal  guidance  law  is  obtained.  This  guidance  law  is  simpler  and 
can  make  the  terminal  trajectory  steep  too.  Therefore,  it  could  be  applied  to  projects  more  easily  and 
requires  less  onboard  computational  resources.  However,  it  consumes  slightly  more  control  energy  than 
the  optimal  guidance  law. 

©  201 1  Elsevier  Masson  SAS.  All  rights  reserved. 


1.  Introduction 

Global  positioning  system  (GPS)-based  guidance  systems  have 
fire  and  forget  feature,  are  weatherproof  and  low  cost  per  kill, 
and  have  therefore  received  much  attention  from  researchers  [13, 
1-4,11],  The  vertical  measurement  error  of  GPS  is  larger  than  its 
horizontal  error  [9],  In  order  to  reduce  its  effect  on  guidance  ac¬ 
curacy,  a  near-vertical  descent  is  required  in  the  terminal  guidance 
phase.  The  terminal  flight-path  angle  is  expected  to  be  close  to, 
but  not  equal,  —90°  because  the  constraints  of  volume  and  cost 
limit  the  control  capability  of  guided  munitions.  The  large  landing 
angle  can  also  increase  the  damage  effect  and  penetration  capacity 
of  guided  munitions.  Previous  guidance  law  research  related  to  an¬ 
gle  constraints  may  be  largely  divided  into  two  categories:  optimal 
[15,10,14,8]  and  not  optimal.  The  latter  is  usually  more  likely  to 
have  some  advantages,  such  as  having  a  simpler  form,  not  requir¬ 
ing  range-to-target  information  [12],  not  requiring  a  time-to-go  tg0 
estimation  [7],  and  being  robust  [6,5].  To  optimize  the  combination 
of  landing  angle,  miss  distance,  and  control  energy  consumption, 
research  is  needed  regarding  the  optimal  guidance  law  with  multi¬ 
constraints. 

In  order  to  design  the  optimal  guidance  law,  motion  equations 
should  first  be  established.  These  equations  often  include  time- 
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varying  parameters.  However,  for  time-varying  systems,  the  opti¬ 
mal  guidance  law  cannot  generally  be  solved  analytically.  There¬ 
fore,  time-varying  parameters  are  not  intended  to  appear  in  the 
design  of  guidance  law.  When  time-varying  parameters  are  in¬ 
evitable,  there  are  two  methods  that  may  solve  the  problem.  One 
method  is  substituting  constants  for  time-varying  parameters.  This 
method  is  simple,  and  is  therefore  employed  by  many  studies; 
however,  the  simplification  critically  decreases  the  effect  of  the  de¬ 
signed  guidance  law.  In  fact,  this  method  cannot  be  applied  when 
the  variation  range  of  the  time-varying  parameters  is  large.  The 
other  method  is  combining  the  optimal  control  theory  and  the 
numerical  value  method.  However,  the  method  has  not  been  ac¬ 
tively  reported  in  the  design  of  optimal  guidance  law  with  multi¬ 
constraints. 

Block  pulse  function  series  are  orthogonal  functions  with  a 
piecewise  constant  value.  They  are  simple  and  easy  operations. 
In  this  paper,  a  block  pulse  function  series  is  applied  in  order  to 
solve  the  problem  described  above.  For  time-varying  control  sys¬ 
tems,  a  design  method  is  proposed  for  the  optimal  guidance  law 
with  multi-constraints  combining  the  numerical  value  method  and 
the  optimal  control  theory. 

2.  Motion  equations 

Consider  the  homing  guidance  geometry  shown  in  Fig.  1.  Here, 
M  and  T  denote  the  guided  projectile  and  target,  respectively,  and 
they  are  regarded  as  particles  whose  movements  are  in  a  vertical 
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Fig.  1.  Homing  guidance  geometry. 


plane.  Take  the  points  M,  MF,  and  r  =  MT  as  the  pole,  polar  axis, 
and  projectile-target  distance,  respectively.  V  is  the  projectile  ve¬ 
locity,  and  0  is  its  flight-path  angle.  GPS-guided  munitions  are  used 
against  stationary  targets;  thus,  the  target  velocity  is  zero.  The  po¬ 
lar  equations  of  motion  in  this  homing  problem  are  given  by 

r  = -V  cos(q  -  0),  and  (1) 

rq=V  sin(q  -  0 )  (2) 

Differentiating  (2)  with  respect  to  time  and  substituting  (1)  into  it, 


as  possible,  and  the  landing  angle  should  be  close  to  90°.  There¬ 
fore,  the  performance  index  can  be  defined  as  follows: 

tf 

J  =  XT(tf)CX(tf )  +^J  R(t)u2(t )  dt  (9) 

to 

where  C  =  [  J  °],  R(t)  =  1,  and  tf  =  to  +  Based  on  the  optimal 
control  theory  [17],  the  optimal  guidance  law  for  the  system  above 

U(t)  =  K(t)X(t )  (10) 

where 

K(t)  =  -R^BT(t)P(t)  (11) 

Here,  K(t)  is  the  guidance  law  coefficient  matrix  that  will  be  ob¬ 
tained  by  the  following  solution  procedure  and  P  satisfies  the 
Riccati  matrix  differential  equation: 

P(t)  =  -P(t)A(t)-AT(t)P(t) 

+  P(t)B(t)R-\t)BT(t)P(t)  -  Q.  (t)  (12) 

where  Q  (t)  =  [°  jj],  and  the  terminal  condition  of  P(t)  is 
P(tf)  =  C  (13) 


Eq.  (12)  can  be  decomposed  into  the  following  two  matrix  differ¬ 
ed)  ential  equations: 


In  general,  GPS-guided  munitions  have  little  or  no  thrust.  There¬ 
fore,  they  cannot  maintain  a  constant  axial  velocity,  and  their  ve¬ 
locity  changes  slowly.  Thus,  V /V  1  and  V /V  =  0  is  assumed. 
Thus,  (3)  can  be  written  as 


W(t)  =  /l(t)  W(t)  -  B(t)R~i  (t)BT(t)Y(t)  (14) 

Y(t)  =  -Q(t)W(t)-AT(t)Y(t)  (15) 

The  terminal  conditions  of  1/1/  (t)  and  Y  (t)  are 


q=-2rrq+rr6  (4) 

Let  Xi(t)  =  q  -  0p,x2(t)  =  Xi  (t)  =  q,  and  u(t)  =  9.  Here,  9F  is  the 
expected  flight-path  angle  in  landing  moment  tf,0F  =9(tf),  and 
u(t)  is  the  guidance  law  that  we  want  to  obtain.  Therefore,  from 


(4),  we  have 

Xi  (t)  =  X2(t)  (5) 

X2(t)  =  — 2^x2(t)+  ~u(t)  (6) 

Let  g(t)  =  — r/r.  Therefore,  from  (6),  we  have 

k2(t)  =  2g(t)X2(t)-g(t)u(t)  (7) 

Based  on  (5)  and  (7),  the  motion  equations  can  be  written  as 

X(t)  =  A(t)X(t)  +  B(t)u(t)  (8) 

where 


x(t)=h(t)i  A(t)-r°  1 1  BityJ  0 1 

L*2(t)J  Lo  2g(t)J 

Initial  conditions:  when  t  =  to,  Xi  (to)  =  q(to)  -  9f,  and  x2(to)  m 

q(to)- 

3.  Optimal  guidance  law  design  on  the  basis  of  block  pulse 
functions 

In  the  guided  process,  the  control  energy  consumption  should 
be  as  small  as  possible,  the  terminal  projectile  should  be  as  steep 


W(tf)  =  I  (16) 

Y(tf)  =  C  (17) 

where  /  is  the  second  order  unit  matrix.  The  solution  for  the  Ric¬ 
cati  equation  is: 

P(t)  =  Y(t)W~\t)  (18) 


Rearrange  (14)  and  (15)  as  follows: 


Z(t)  =  F(t)Z(t) 
where 


Z(t)  = 


F(t)  = 


1/1/  (t)  1 

.  Y(0  \ 
A  (t) 
.-0.(0 


— BCtjU-1  (t)HT(t)  1 

-AT(t)  J 


(19) 


(20) 

(21) 


Integrating  (19)  from  tf  to  t,  we  can  derive  the  following: 


Z(t)-Z(tf)  =  j  F(r)Z(r)dr  (22) 

tf 

Expanding  (22)  into  block  pulse  series,  we  have 


[ z !  Z2  ■■■  Z(tf)  •••  Z(t/)]#(t) 

=  -[F,Zx  F2Z2  ■■■  FmZm]hHTd>(t)  (23) 

where 
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^(0  =  [0l(t)  <h  (!)  </>m(t)]T 


(24) 

(25) 


Here,  4>i(t)  is  a  block  pulse  function,  i  =  1, 2, . . . , m;  m  is  the  num¬ 
ber  of  block  pulse  functions,  and  h  is  the  block  pulse  width.  Fk 
is  the  mean  value  of  F(t)  in  the  interval  [(k  —  1  )T /m,  kT /m].  The 
definition  and  some  elementary  properties  of  block  pulse  functions 
can  be  found  in  Appendix  A.  From  (23),  we  have 


Zr  -  Z(tf )  =  -~h(F  1Z1  +  2F2Z2  +  •  •  •  +  2FmZm)  (26) 

Z2  -  Z(tf )  =  -^h(F2Z2  +  2F3Z3  +  •  •  •  +  2FmZm)  (27) 


Zm_i  -Z(t/)  =  --h(Fm_1Zm_1+2FmZm)  (28) 

Zm-Z(t/)S-i/l(FmZm)  (29) 


Solving  Eqs.  (26)-(29),  we  obtain  the  following  recursive  algo¬ 
rithm: 


m=(j+lhFm)  Z(tj) 


(30) 


1, 


From  (20),  (30),  and  (31),  we  can 
1, 2 . m.  Eq.  (11)  can  be  written 


K  (t)  =  [  if  i  K2 


/--hF^JZfc,  k=2,3, . . .  ,m 

(31) 

obtain  Wj<(t)  and  Yk(t),  k  = 


Fig.  3.  k2  history  vs.  tg0. 


1000- 

1  °' 

(32) 

5000  10000  15000  20000 


where 


Ki  =  [khi  k2,i  ]T ,  i=  1,2, . . .  ,m  (33) 

Rk  is  the  mean  value  of  R(t)  in  the  interval  [(k  —  l)T/m,kT/m], 
and  B[  is  the  mean  value  of  B(t )  in  the  interval  [(/-  1)T /m,  IT /m]. 
Based  on  the  orthogonality,  (32)  can  be  written  as 

K(t)  =  -J2Rk'BTkYkW~'Mt)  (34) 

k=l 


4.  Simulation  results 

In  a  ground  coordinate  system,  let  the  muzzle  velocity  of  the 
guided  projectile  be  800  m/s,  and  its  initial  flight-path  angle  do 
is  30°.  The  muzzle  coordinates  are  (0, 0),  and  the  target  coordi¬ 
nates  are  (20000,0).  The  guidance  will  begin  when  the  projectile 
reaches  the  highest  point.  Figs.  2  and  3  show  the  histories  of  ki 
and  k2  versus  tg0,  respectively,  where  tg0  —  tf  -  t.  Ballistic  sim¬ 
ulations  were  conducted  using  three  lands  of  guidance  laws,  and 


Fig.  4.  Trajectories  using  three  different  guidance  laws,  respectively. 

the  results  were  compared.  The  first  guidance  law  is  from  the  pro¬ 
posed  method,  called  the  optimal  guidance  law.  The  second  is  the 
suboptimal  guidance  law,  which  is  obtained  by  taking  the  means 
of  ki  and  k2  as  guidance  law  coefficients.  The  third  is  the  propor¬ 
tional  navigation  law. 

Fig.  4  compares  the  three  trajectories  achieved  with  the  three 
guidance  laws  above.  Lines  1,  2,  and  3  are  computed  using  the 
optimal  guidance  law,  suboptimal  guidance  law  and  proportion 
navigation  law,  respectively.  The  same  is  indicated  in  the  follow¬ 
ing.  Fig.  2  shows  that  k i  increases  as  time-to-go  decreases,  and 
Fig.  3  shows  that  fc2  decreases  as  time-to-go  decreases.  As  men¬ 
tioned  above,  the  coefficients  of  the  suboptimal  guidance  law  are 
the  means  of  ki  and  k2.  In  the  first  phase,  the  mean  of  ki  is  2-3 
times  the  original  value,  and  in  the  final  phase,  it  is  1/2- 1/5  of  the 
original  value.  Fig.  4  shows  that  the  trajectory  2  is  a  good  way  to 
improve  positioning  accuracy  because,  for  terminal  guidance,  the 
trajectory  is  steep,  though  needing  more  control  energy. 


204 


L  Dou,J.  Dou  /  Aerospace  Science  and  Technology  23  (2012)  201-205 


Table  1 

Comparison  of  three  guidance  laws. 

Guidance  law  Landing  angle  Control  energy  consumption 

(degree)  (substituted  with  the  absolute 

value  integral  of  over  loading) 

Optimal  guidance  law  84.85  26.27 

Suboptimal  guidance  law  80.74  30.10 

Proportional  navigation  38.94  12.47 


Miss  distance  (m) 


f/s 


paper  attempts  to  solve  this  problem  by  combining  the  block  pulse 
functions  with  the  optimal  control  theory.  Expanding  the  time- 
varying  systems  into  the  block  pulse  series,  the  system  can  be 
translated  into  a  series  of  time-invariant  systems,  which  can  be 
solved  with  the  optimal  control  theory.  The  results  show  that,  us¬ 
ing  the  proposed  optimal  guidance  law,  the  landing  angle  is  three 
times  as  large  as  what  is  obtained  when  using  the  proportion  nav¬ 
igation  law.  The  steep  trajectory  can  increase  strike  accuracy,  not 
only  because  it  increases  the  GPS-based  guidance  accuracy,  but 
also  because  it  decreases  the  error  caused  by  time  step.  The  sub- 
optimal  guidance  law  has  a  very  simple  form  and  can  save  signif¬ 
icantly  computation  resources;  however  it  consumes  more  control 
energy. 


Fig.  5.  Flight-path  angle  history. 


Acknowledgements 


t/s 


Fig.  6.  Absolute  value  integral  of  overload. 


This  work  was  supported  by  National  Natural  Science  Foun¬ 
dation  of  China  (Grant  No.  60904085),  New  Teachers’  Fund  for 
Doctor  Stations  of  Ministry  of  Education  of  China  (Grant  No. 
200802881012),  “Excellent  Talent  Project  Zijin  Star"  Foundation  of 
Nanjing  University  of  Science  and  Technology,  Foundation  of  Na¬ 
tional  Defence  Key  Laboratory  of  Ballistics. 

Appendix  A 


Definition  of  block  pulse  functions  and  their  some  elementary 
properties  [16]. 

A  block  pulse  function  set  0fc(t)  (k  =  1, 2, . . . ,  m)  can  be  defined 
in  the  interval  [0,  T]  as 


J  1,  (k-1)T/m^t  <kT/m 
1 0,  otherwise 


(35) 


The  landing  angles,  control  energy,  and  miss  distances  using  the 
three  different  guidance  laws  are  presented  in  Table  1.  Here,  the 
purpose  is  to  compare  control  energy  consumption  focusing  purely 
on  the  relative  comparisons.  Therefore,  the  absolute  value  integral 
of  overload  is  substituted  for  the  control  energy  consumption.  The 
overload  is  defined  as  the  ratio  N/G,  where  G  is  the  gravity  acting 
on  the  guided  projectile  and  N  is  the  sum  of  all  external  forces, 
except  gravity,  acting  on  the  guided  projectile.  The  miss  distance  is 
defined  as  the  straight-line  distance  between  the  landing  point  of 
the  guided  projectile  and  the  target.  Table  1  shows  that  the  con¬ 
trol  energy  consumption  using  the  suboptimal  guidance  law  is  15% 
higher  than  what  is  needed  when  using  the  optimal  guidance  law. 
Figs.  4-6  show  that,  when  using  the  proportion  navigation  law,  the 
control  energy  consumption  is  smallest;  however,  the  landing  an¬ 
gle  is  flattest  and  only  38.94°.  Therefore,  the  proportion  navigation 
law  cannot  satisfy  the  trajectory  demand  of  GPS-guided  munitions. 

5.  Conclusions 

In  order  to  optimize  the  combination  of  landing  angle,  miss 
distance,  and  control  energy  consumption,  research  regarding  the 
optimal  guidance  law  with  multi-constraints  is  necessary.  However, 
the  optimal  guidance  law  cannot  generally  be  solved  analytically 
when  the  motion  equations  include  time-varying  parameters.  This 


The  block  pulse  functions,  (</>o(t),  0i(t), . . . ,  0m(t)},  are  dis¬ 
joined  from  one  another  in  [0,  T] : 

.  \  M),  k  =  j 

0fc(O0j(t)  =  |  o  (36) 

The  block  pulse  functions  are  orthogonal  to  one  another  in  [0,  T] : 
}  fT/m,  k  =  j 

j  &(0*j(t)dt=  jQ/  k^J.  (37) 

o 

An  arbitrary  integrable  function  /(t)  in  [0,  T]  can  be  approxi¬ 
mately  expanded  into  the  block  pulse  series  as 

=  (38) 

k=  1 
Let 

*(t)  =  [fc(t)  02(f)  •••  0m(f)]T  (39) 

The  integral  of  the  block  pulse  series  can  be  written  as  follows: 

f  0(t)dt=-H0(t),  0  <t<T  (40) 

J  rn 
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where  H  is  the  m  x  m  matrix: 
12  2  •••  2" 

0  1  2  •••  2 

0  0  1  •••  2 

.0  0  0  •••  1. 
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